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I. INTRODUCTION 



Accurate solutions for discrete states of few-nucleon systems interacting via a realistic 
potential have been obtained with various sophisticated methods pj. The most significant 
merit of these methods is that the interaction employed is tested strictly in comparison with 
experiment and thus the nuclear structure can be understood without ad hoc assumptions, 
and some important ingredients such as a three-body force are revealed. In contrast to the 
bound-state problem, a microscopic description of reactions using a realistic potential has 
been hampered by the difficulties related to continuum states as well as scattering boundary 
conditions. Some progress has recently been made towards the application of an ab initio 
approach to the problem of scattering and resonances @, Q. The study in this direction 
should be further pursued as it may provide more detailed information on the characteristics 
of the interaction as functions of partial waves and incident energies and on the dynamics 
of the participating nuclei. 

Because of their apparently different natures, continuum and discrete states are obtained 
separately using different methods. Since, however, both of them are solutions of the same 
Schrodinger equation, it would be nice if they could be obtained on an equal footing, namely 
if the correlation of the constituent particles could be included in the scattering problems 
as well as in the bound-state problems. Several methods have in fact been proposed us- 
ing square- integrable (£ 2 ) basis functions for continuum problems in atomic and nuclear 

physics BSS^HB- 

The purpose of this article is to show a simple approach to the scattering problem using 
the technique for discrete states, particularly developed for a basis expansion method. The 
present approach has some similarity to that of Ref. |6|, where the stochastic variational 



method |10l . Illl ] is used in conjunction with stabilization ideas to calculate the low energy 
phase-shifts for positronium-atom scattering. The basis set for describing the configuration 
space of the interaction region is spanned by the stochastic variational method. In Ref. 0], 
the phase shift is extracted by fitting the wave function in the scattering region to sin(fcr + (5). 
We instead use the Green's function approach to describe the wave function in the scattering 
region. 

The basic quantity in the scattering of nuclei is the overlap integral between the product 
of the two internal wave functions of the nuclei and the scattering wave function of the 
composite system. This overlap integral is called a spectroscopic amplitude (SA) in this 
article. The phase shift is determined from the asymptotic behavior of the SA. This type 
of overlap integral was studied long time ago for extracting spectroscopic information from 



single-particle transfer reactions 12|, [l3|. In Sec. [Ill we discuss the equation of motion which 
the SA satisfies, and we derive a formula to calculate the phase shift. The accuracy of the 
present approach is tested in Sees. IIII Al and IHI Bl Section IIII CI presents an application to 
a + n scattering using a realistic nucleon-nucleon potential model. A brief summary is given 
in Sec. [TV] The calculation of the SA is explained in the Appendix. 



II. FORMALISM 

Let c stand for a channel including a pair of nuclei denoted cii and cti whose normalized 
wave functions, V^it^i) an d ipi 2 (a 2 ), are assumed to be given. The wave function of the 
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total system with energy E generally takes the form 



c 7 

The first term on the right side stands for configurations for the elastic channel as well as 
other channels, whereas the second term represents the configurations that are needed to 
take account of the effects of resonances, distorted states, etc., not included in the first term. 
The term S 7 jjv/ is thus assumed to represent the configurations in which all the nucleons 
are confined in the interaction region. It is assumed that \& jm satisfies the antisymmetry 
requirement for the exchange of the nucleons. 
A key in the present method is the SA 

y(r) = ($ cJM {r)\*j M ), (2) 
where $ c jM( r ) is a test function for the channel c 

**w(r) = [[^>i)^> 2 )]/YHf c )] JM 6{r l~ r) . (3) 

The coordinate r c denotes the relative distance vector of the two nuclei. The angular 
momenta l\ and I2 of the two nuclei are coupled to the channel spin J, which is then 
coupled with the orbital angular momentum i of the relative motion to the total angular 
momentum JM. In the test function Q, the two nuclei are pinned down at the relative 
distance specified by r. 

The phase shift for the scattering of nuclei a% and «2 is calculated from the asymptotic 
behavior of y(r). When ^ jm is obtained in a bound-state approximation, the SA calculated 
from it has ill behavior at large distances, and hence the asymptotics of y(r) is usually not 
good enough to enable one to calculate the phase shift accurately. To resolve this problem, 
we derive an equation of motion which y{r) calculated using the exact ^ jm should satisfy. 
For this purpose we start from the equation 

(QcJM(r)\H\*j M ) = E(<S> cJM (r)\y JM ). (4) 

The Hamiltonian of the system can be decomposed into 

H = H ai + H a2 + T C + V e , (5) 

where H ai and H a2 are, respectively, the internal Hamiltonians of nuclei ol\ and ct2, 

h 2 d 2 



is the kinetic energy of the relative motion between them, and 

v c = £ % (7) 



is the interaction acting between the two nuclei. If the Hamiltonian contains a three-body 
force, the potential V c acting between the nuclei should include the following terms of the 
three-body potentials, E( ieai) <(i<fce«2) v m+ Y,( l<j&ai )<(k^ 2 )^k- 
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Substituting Eq. ((SJ) into Eq. (J4]) and introducing a local potential U c (r) which acts be- 
tween nuclei ot\ and «2 makes it possible to transform Eq. (j3J) to the following inhomogeneous 
differential equation for y(r) (the subscript c in C/ c and /i c is suppressed hereafter) 



d 2 2d 
dr 2 r dr 



£{£ + 1) 2 fx 



-§U{r) + k 2 



(8) 



where k—y2/i(E — E ai — E a . 2 )/Ti 2 with = {ip IxMi \H a .\il)i. M .) is the wave number for the 
relative motion, and the functions z(r) and w{r) are, respectively, defined by 



and 



z(t) = {$ b jm{t)\V c -U\Vj M ), 



w{r) = ($ cJA /(r) | H ai - E ai + H a2 - E a2 \ <S> JM ). 



(9) 



(10) 



The function w(r) vanishes if ^m, are the eigenfunctions of H a ., that is, H ai il) IiM .—E a .ip I . Mi . 

Equation (jSJ) is apparently equivalent to Eq. PJ), which is the Schrodinger equation pro- 
jected to the space spanned by the test function $ c jj\/(r). The equivalence does hold for an 
arbitrary choice of U(r). Let v(r) and h(r) denote, respectively, the regular and irregular 
solutions of the homogeneous equation with z(r) + w(r) being set to zero in Eq. dSJ). They 
are defined to satisfy the Wronskian relation, W(v, h)(r) = v(r)h'(r) — v'(r)h(r)=l/(kr 2 ). 
A general solution of Eq. (jSJ) that is regular at r=0 and has the asymptotic behavior appro- 
priate to the scattering solution reads 



2/i 

H .7 



y(r) = Xv(r) + 7^ / G(r,r')[z(r') +w(r')]r' dr' 



(11) 



where A is a constant. Here the Green's function G which is a solution of the following 
equation 



d 2 2d 
dr 2 r dr 



+ 1 



pl/(r) + k 2 



1 



G( r ,r') = —5{r-r') 



is given by |1 



G(r,r') 



kv{r)h{r') r < r' 
kh(r)v(r') r > r' . 

By expressing the integral in Eq. (TTTj) as 



G(r,r')[z(r') +w(r')]r' dr' = k[p(r)h(r) + q(r)v (r)] 



with 



p(r) 



v(r')[z(r') + w(r')]r' 2 dr' 



q(r) 



h{r')[z{r') + w(r')]r' 2 dr' 



(12) 



(13) 



(14) 



(15) 



the SA (TIT]) takes the form 



A + 



2fjtk 



q(r) 



v(r) H — ^—p(r)h(r). 

h 



(16) 



4 



Taking the asymptotics of this y(r) determines the phase shift 5e as 



tanSe = tan^ -* — -^-p(oo), (17) 

where 8g is the phase shift corresponding to the potential scattering by U(r). Let y(r) of 
Eq. (USD (or Eq. (JT])) be called SAGF (SA solved with the Green's function). 

Equation ( fTTj) shows that we can obtain the phase shift accurately if A and p(oo) are 
known to high accuracy. The value of p(oo) consists of two terms, one involving the function 
z(r') and the other involving the function w(r'). As was already mentioned, w(r') vanishes 
if ipiiMi are the eigenfunctions of the internal Hamiltonian H ai . As we will see later, even 
though ipj i M i are not completely identical with the eigenfunctions, the magnitude of w(r') 
turns out to be much smaller than that of z(r'). The function z(r') consists of the sum 
of various pieces of V c such as central, tensor and spin-orbit forces, and thus tan 5e can be 
decomposed into the contributions of those terms. 

In a practical calculation of a phase shift, we have only an approximate solution for *&jm, 
and because of this y(r) as well as z(r) and w(r) are all approximately evaluated. The value 
of A is determined by comparing the two SAs, Eqs. (j2I) and ([TBI) , as will be discussed later. 
Evaluating p(oo) to a good approximation depends on a choice of U(r). We assume that 
U is chosen in such a way that V c approaches U for large r, namely V c — U is an operator 
that is non-zero only in the interaction region. If U is chosen to satisfy this condition, z(r) 
can be accurate even though ^ jm does not have a correct tail in the region where V c — U 
is negligibly small. The function w(r) is accurate as well in the interaction region. Then 
we may assume that both p(r) and q(r) can be evaluated fairly accurately provided U is 
suitably chosen. 

The functions z(r) and w(r) (and y(r) of Eq. fl2])) are calculated using a code for bound- 
state calculations in so far as the Dirac 5-function in the test function ([3]) is approximated 

as 

^-^-E/.w.w, ( 18 ) 

r c r z — ' 

where {f u } is an C? 'pseudo-complete' set with (/„ | f u >) = S vu i [15j. The expansion like 
Eq. ([TBI) is used also in Ref. J3] in the basis of harmonic-oscillator functions. Because a more 
precise evaluation of y(r) is in general desirable, we show in the Appendix an analytical 



method to calculate the SA pj) for the correlated Gaussian basis functions [UJ, [16| which 
are employed in this article. The accuracy of the series expansion ( {TBI) in the calculation 
of functions of the type of z(r) was discussed in Ref. I3| by comparing it to the exact 



calculation for the a + 1 system where the two fragments are assumed as the lowest shell- 
model states with a common oscillator parameter. 

An equation similar to Eq. (J5J) was proposed to improve the overlap integrals which 



appear in nucleon transfer reactions or virtual nucleon decay [18|, [19j and a decay |2C 
To our knowledge, no one has yet used Eq. (jSj) to obtain the phase shift in a microscopic 
calculation. 

In a microscopic reaction theory such as the resonating group method (RGM) [2lj], ^cjM 
of Eq. flTJ is expressed as 

POD 

VcJM= / u e (r)A$ cJM (r)r 2 dr, (19) 
Jo 
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where A is the internucleus antisymmetrizer, and the phase shift for the elastic scattering is 
calculated from an integro-differential equation for ue(r) [22]. We stress that y(r) is used in 
the present approach instead of ue(r). The function y(r) is subject to the simple differential 
equation, and it is determined uniquely regardless of whether or not Pauli-forbidden state 
exist, which is in contrast to the case of ui{r). 

The accuracy of the phase-shift calculation crucially depends on how accurately A is 
determined. We tested two ways to determine A. The first is to fit the SAGF (1161) which is 
a function of A to the SA (T5]) with the least squares method in the interval [r , n] where the 
SA (T5]) is expected to be accurately obtained: 

minimize over A : V" [y SAGF (?"i) - Z/ SA (^)] 2 - (20) 

i (r <r 4 <ri) 

We found that A determined in this way remains virtually unchanged within moderate 
choices of the interval. The second is to calculate the Wronskian W(y,h)(r) using both 
y SA (r) and y SAGF (r). The latter reads W(y SAGF , h)(r) = (\/kr 2 ) + (2[i/h 2 r 2 )q(r). Equating 
the two Wronskians leads to the following expression for A: 

A(r) = kr 2 W{y SA , h){r) - ^g(r) = kW(ry SA , rh){r) - ^q(r), (21) 

n h 

which usually becomes r-dependent because y SA (r) is only approximately equal to y SAGF (r). 
The least squares fitting to this A(r) in the interval [r ,ri], 

minimize over A : [A(?"i) — A] 2 , (22) 

i (ro<ri<ri) 

yields again a stable A value, which is in good agreement with that determined by the first 
method. We use the first one mostly in what follows because it requires no differentiation 
of y(r). Interestingly, unlike the i?-matrix theory j23|, our phase-shift calculation requires 
no channel radius. 

As shown in Ref. [3], the construction of the Green's function is easy even for coupled- 
channel problems if U is local. It is thus noted that the present approach can be straight- 
forwardly extended to the scattering including coupled-channels. 



III. EXAMPLES 

A. 3 5i n + p scattering 

A first example to test the present approach is the 3 Si n + p scattering phase shift 
calculated with the Minnesota potential [24|. This is just a potential scattering of the two 
particles, and a numerically exact phase shift can easily be obtained. Diagonalizing the n+p 
Hamiltonian with the Minnesota potential v(r) in appropriate C 2 basis functions produces, 
besides the deuteron ground state, continuum discretized states and corresponding energies 
E. The phase shifts are calculated using these states. The C 2 basis functions used for the 
S wave are Gaussians, exp(— \(3r 2 ) with different falloff parameters f3, where (3 is real or 



complex with Re/3 >0 25|. In this potential problem, y SA (r) for the discretized energy E 
takes the form 



K 

VT\r 



SA (^-J2d(E)e~ ¥ir \ (23) 



i=l 
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TABLE I: Comparison of the 3 Si phase shifts, given in degrees, of n+p scattering between Numerov 
and SAGF methods. The A value is determined using Eq. (J20J) (Method 1) or ([22} (Method 2) 
with the use of different intervals [ro,ri] (fm). The Minnesota potential 24| is used. 
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and the function z(r) reduces to v(r)y SA (r), and w(r) vanishes. The potential U(r) is set 
to zero. 

Table fl] compares the phase shifts calculated using Eq. ( TT7T) with those obtained with the 
Numerov method, which is virtually exact. The comparison shows that the present method 
produces very stable phase shifts which are rather insensitive to the choice of the method of 
determining A as well as the interval used for the minimization of the error. To generate the 
discretized states for different energies, we repeated the calculation by changing K and the 
set of (Pi, p 2 , ■ ■ ■ ,Pk)- Figure [1] displays the phase shifts obtained in this way. The phase 
shifts obtained with our method almost perfectly agree with those of the Numerov method 
in a wide range of incident energies. 
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FIG. 1: Comparison of the 3 5*1 phase shifts of n+p scattering between SAGF and Numerov meth- 
ods. The Minnesota potential [24| is used. 



Numerov 




i 



7 



B. a + n scattering with an effective nucleon-nucleon potential 



In this and the following subsections we calculate the S- and P-wave phase shifts of 
the a + n scattering in a microscopic approach. The nucleon-nucleon interaction employed 
in this subsection is the Minnesota potential (24| which consists only of the central and 
spin-orbit potentials. The u parameter of the central potential is set equal to 0.98, and the 
spin-orbit potential adopted is — 591.1 e _3r £ ■ a in MeV (set IV of Reichstein and Tang). 
This potential is known to reproduce the empirical a + n phase shifts of Refs. j^jj , 27| . The 



phase-shift difference between the 3/2" and l/2~ states is particularly well reproduced by 
this spin-orbit potential, and in this sense the Minnesota potential can be regarded as an 
effective potential. 

Only a single a + n channel is included in the phase-shift calculation. The wave function 
of the a particle, ipo(a), is obtained by diagonalizing the Hamiltonian H a in a basis of a 
number of Gaussians. The binding energy with the Coulomb potential being included is 
29.90 MeV and the root-mean-square matter radius is 1.41 fm. 

The relative motion function U((r) in Eq. (I19p is taken as a combination of Gaussians: 

K 

U£ (r) = J2Cir e e-^ r \ (24) 

i=l 

We have calculated the phase shifts using two different methods, the microscopic i?-matrix 
theory (28| and the present method. The accuracy of the /2-matrix theory is well tested, 
and it is considered to produce virtually exact results to which the phase shifts of SAGF 
are to be compared. In the SAGF calculation the parameters bi = with real $ are 

chosen to form a geometric progression to cover < bi ^ 10 fm. The number K of the 
basis functions is about 10-15. The Gaussian basis used in Eq. ( 1241 can also be employed to 
construct the pseudo-complete set in Eq. f|T8|) . We obtained all the matrix elements needed 



in SAGF using the method developed in Ref. [16|. See also the Appendix for the calculation 
of the SA. 

Figure [2] compares the SAGF phase shifts with the i?-matrix phase shifts. The value of A 
is determined by the first method using the interval [1,5] (fm). The agreement is excellent. 
Both of z(r) and w(r) are included in the SAGF calculation. We switched off the function 
w(r) and found that the phase-shift change is negligible: The largest change of a few percent 
occurs in the resonance region at around 1 MeV of the 3/2~ phase shifts. Except for this case 
the phase-shift change is smaller by one order of magnitude. Thus we may safely neglect 
the contribution of w(r). 

It is interesting to see how much the SAGF changes from the SA. Figure [3] compares the 
two SAs for the 3/2~ state at three different energies. In each case, the SA curve agrees 
very well with the corresponding SAGF in the region of r < 5fm, which indicates that the 
expansion (|2~41) is good enough to describe to a good approximation the a-n relative motion 
function in the interaction region. The peak position and the amplitude of the SA curve 
begin to deviate from those of the SAGF curve for r > 5 fm. The deviation becomes larger 
as the energy increases. 
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C. a + n scattering with realistic nucleon-nucleon potentials 



In this subsection we take the same model as in the previous subsection but employ the 
realistic nucleon-nucleon potentials of AV8 type, AV8' 2^] and G3RS [30]. Both of them 
contain central (V c ), tensor (Vt) and spin-orbit (Vb) terms. The L 2 and (L ■ S) 2 terms of 
the G3RS potential are ignored. The binding energy of the a particle is 25.09 for AV8' and 
25.29 MeV for G3RS 0. In what follows we use slightly truncated wave functions for the 
a particle to save computer time. The truncation is done by ignoring the small component 
with total orbital and spin angular momenta L = 1, S = 1 the magnitude of which is of 
the order of 0.3% or by optimizing the wave function in a smaller basis set whose size is 
approximately one third of the converged solution. The change of the phase shift due to the 
truncation is estimated to be at most a few percent. 

Figure H] displays the S- and P-wave phase shifts obtained in an a + n single-channel 
calculation. The difference due to the nucleon-nucleon potential is not very large. The 
agreement between theory and experiment is quite reasonable for the l/2 + phase shifts. Our 
phase shifts agree fairly well with those of the Quantum Monte Carlo calculation obtained 
using only the two-body potential of AV18 [2|. In a sharp contrast to the 5-wave phase shift, 
the calculated P-wave phase shifts are by far smaller than experiment, and considerably 
smaller than those obtained in Ref. 0]. We have coupled an inelastic channel of 4 He(0j) +n 
with the elastic channel, but the effect of coupling is negligible below E =15 MeV. 

In Fig. [5] we show the contribution of the potential components to the S'-wave phase 
shift, tan5 - Generally speaking, the two potentials give similar results. Because the spin- 
orbit force produces a negligible contribution, we do not show its contribution in the figure. 
The phase shift gets largest contributions from the central and tensor forces. For the AV8' 
potential, the tensor contribution is larger than the central contribution, whereas, for the 
G3RS potential, the central force gives a larger contribution. This is quite consistent with 
the relative importance of the tensor and central forces found in the binding energies of 
A = 3, 4 nuclei fTF 



1 80 i 1 1 1 1 1 r 




E [MeV] 



FIG. 2: Comparison of the phase shifts of a+n scattering between the microscopic i?-matrix and 
SAGF methods. Solid, dashed and dotted lines are the phase shifts calculated with the i?-matrix 
theory, while symbols are those by SAGF. The Minnesota potential [3] is used. 
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FIG. 3: Comparison between the SA and the SAGF for the a+n scattering (3/2 ) at three different 
energies. The Minnesota potential 2J] is used. 



The experiment shows a 3/2~ sharp resonance around E = 0.9 MeV and a l/2~ broad 
resonance. The single-channel calculation apparently misses some important configurations 
which are relevant to the resonances, for example the inelastic channels including the excited 
states of 4 He with J 71 ^ + , as indicated in Ref. {3j. Other configurations which could be more 
important than these inelastic channels include the channels of different cluster partitions, 
3 H+d and 3 He+2n, where 2n denotes a di-neutron cluster. The threshold energies of these 
channels are lower than those of the inelastic channels. The clusters 3 H, d, and 3 He may 
not necessarily be in their ground states but can be in their pseudo-excited states. This 
is plausible because the strong tensor component included in the realistic potential brings 
about a mixing of different orbital and spin angular momenta of the participating nuclei. 

When we include the different cluster partitions of 3 H+d and 3 He+2n, we have to note 
the old, unresolved problem that the threshold energies of different channels are not usually 
reproduced well. The difference of the threshold energies between a+n and 3 H+d is 17.6 MeV 
experimentally, but the theoretical value turns out to be about 15 MeV for both the AV8' 
and G3RS potentials. To quantify the effects of the 3 H+d and 3 He+2n channels properly, 
the inconsistency in the threshold energy should be settled beforehand. The inclusion of 
three-body forces is important in this respect. In the calculation of Ref. 0] including the 
three-body force of the Illinois-2 model, the threshold problem does not appear, and some 
of the distorted configurations are certainly included implicitly. However, it is not very 
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FIG. 4: a+ra scattering phase shifts calculated using the AV8' and G3RS potentials. Solid, dashed 
and dotted lines are the empirical phase shifts 26|, |27|. Circle, diamond and triangle symbols 
denote the calculated phase shifts for l/2 + ,3/2~ and 1/2 - , respectively. 



clear whether the improvement in the phase shifts obtained with the use of the three-body 
force indicates simply the predominant importance of the 3 H+<i or 3 He+2n channel or the 
necessity of including further configurations other than the elastic channels of a + n and 
3 H+<i. Clarifying this point will certainly be very important, but it is beyond the present 
work because such calculations require an expensive computation. 

We here remark on similarity in the phase shifts between a + n and 3 He+p scatterings. 
The most attractive phase shifts of the low-energy 3 He+p scattering occur in the P-wave 
channel with 1 = 1 and J = 2, showing a resonance behavior. This is consistent with the 
fact that the ground state of 4 Li has J n = 2~ . The S"-wave phase shift with 1 = 1 and J = 1 
shows a repulsive behavior similarly to the a + n case. Two of the present authors (K. A. and 
Y. S.) and S. Aoyama have recently calculated the 3 He+p phase shifts in a microscopic R- 
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FIG. 5: Contributions of the components of the nucleon-nucleon potential to tan<5o for S-w&ve 
a+n scattering phase shifts. The contribution of the spin-orbit force is negligibly small. 
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matrix method using the G3RS potential [31]. The result of the single- channel calculation 
is very similar to the a + n case. The S-w&ve phase shifts are accounted for fairly well 
in the 3 He+p single-channel calculation, but the calculated P-wave phase shifts turn out 
to be too small compared to the empirical ones. A calculation of including the channels 
of d(l + , + ) + 2p(0 + ) has significantly improved the discrepancy in the phase shifts and a 
further inclusion of the inelastic channels of 3 He+p has reproduced the empirical phase shifts 
reasonably well. We understand this as follows. In the P-wave scatterings of both a + n and 
3 He+p, the nucleon can penetrate close to the nucleus and the nucleon-nucleus interaction, 
particularly the tensor force distorts the nuclear state. Thus the single-channel assumption 
becomes rather poor. On the other hand, in the 5*-wave case the nucleon cannot come close 
to the nucleus because the two identical nucleons occupying the S-orbit repel the incoming 
nucleon due to the Pauli principle, and thus the nucleus receives little distortion. 

IV. SUMMARY 

We have formulated a method to calculate the scattering phase shift using continuum- 
discretized states. The spectroscopic amplitude calculated from the discretized state is 
compared to that obtained with the Green's function, from which the phase shift can be 
determined. The method has been tested successfully in the cases of n+p and a+n scattering 
where two nucleons are assumed to interact via an effective nucleon-nucleon potential. The 
method has the advantage that it can be used with ease to scattering problems of coupled- 
channels. 

Using a sophisticated wave function for the a particle, we have applied the present method 
to a + n scattering described by realistic nucleon-nucleon potentials. The l/2 + phase shifts 
for the S-w&ve showing a repulsive behavior of the underlying an interaction are in reason- 
able agreement with the empirical phase shifts. The missing attraction needed to reproduce 
the data more perfectly is attributed to the effect of three-body forces We have analyzed 
how much the S'-wave phase shifts are affected by the terms of the nucleon-nucleon poten- 
tial. Both the tensor and central forces contribute significantly but the spin-orbit force has 
a negligible contribution. 

The P-wave phase shifts for the 3/2~ and l/2~ states are too small with respect to 
experiment. This unexpected result indicates that a good reproduction of the a + n phase 
shifts attained using an effective interaction such as the Minnesota potential is an ad hoc 
description and may not be very realistic. Based on the comparison with the result of 
Ref . |2| , the discrepancy in the P-wave phase shifts could be reduced by taking into account 
some distorted configurations which are especially important to form the 3/2~ resonance as 
well as the three-body force. A careful study on the effects of the distorted configurations 
will be important to pin down the contribution of the three-body force. 

We would like to thank R. G. Lovas for his careful reading of the manuscript and useful 
comments. This work has been performed as a part of the Bilateral Joint Research Project 
between the JSPS (Japan) and the FNRS (Belgium). W. H. is supported by a Grant-in 
Aid for Scientific Research for Young Scientists (No. 19-3978) as a JSPS Research Fellow for 
Young Scientists. 
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APPENDIX: CALCULATION OF A SPECTROSCOPIC AMPLITUDE 



In this appendix we show a method of calculating y(r) for the basis functions used in 
this article. Because the matrix element for the spin-isospin parts can be calculated in a 
standard technique, we focus on only the spatial part of the SA. The SA plays an important 
role in discussing the spectroscopic properties of nuclear structure. See Ref. [32J for the case 
of 4 He. 

Let (xi,X2, ■ ■ ■ ,xn^ 2 ) denote the internal coordinates of nuclei ot\ and a 2 . Here N is 
the number of nucleons which make up the two nuclei. The set of N — 2 coordinates is 
simply represented by an iV — 2 column vector x whose ith element is X{. For the sake 
of convenience, let x^-i, instead of r c used in the main text, denote the relative distance 
vector between the two nuclei. Let X stand for aniV-1 column vector whose ith element 
is X{. The first N — 2 elements of X are exactly the same as the elements of x, and the last 
element of X is a;jv-i- 

The spatial part of the test function is assumed to take the form 



®LM L (r) 



F{L A L A )L'{U^U^A ,x) Y e (x N _x) 



(Al) 



where the function F is the correlated Gaussian basis state whose angular part is specified 
by two global vectors lfl 



F(l 3 l 4 )L'M' l (u3,U4,A',x) = exp (-^xA'x^j [yL 3 (u3x)y LA (u A x)] L , M , L , (A2) 

with ye m (r) = r e Y£ m (r). Here A' is a positive-definite, symmetric matrix of dimension 
N — 2, and u% and U\ are iV — 2-dimensional column vectors. To simplify the expressions, 
we use the notation xA'x = J2fj=*i A'^Xi ■ Xj and u^x = Yl!i=i u 3i x i- The square bracket 
in Eq. (1A1I) indicates the angular momentum coupling. It will be convenient to rewrite 
Eq. (lAip as follows: 



$r,w, (r) 



J2(L'M' L £m\LM L ) I ' Y tm (r) F {LzU)l , m , l {u^u^ A',x)5{x n ^ - r)dr. (A3) 



Similarly, the spatial part for the basis function of the composite system, a\+a 2 , is 
constructed from the function 

F(l 1 l 2 )lm l ( u 1i M 2, A, X) = exp (^-^XAX^ [yL 1 (u 1 X)y L2 (u 2 X)] LML . (A4) 

Here A is a positive-definite, symmetric matrix of dimension N — l, and U\ and u 2 are column 
vectors of dimension N—l. The RGM wave function ([19]) with Eq. (12^1) can be expressed 
in terms of the above function (IA4|) . That is, the matrix A becomes block-diagonal and 
contains three blocks, each of which has the dimension of N ai — 1, N a2 — 1,1, respectively. 
Here N ai is the mass number of the nucleus «j. The A/v-ijv-i element corresponds to 
Pi. In addition, one of the global vectors, say, u 2 X must be equal to £Cjv_i and the other 
global vector, U\X, should not contain ajjv-i- This requirement is met simply by choosing 
u 2l — (i — 1, . . . , iV — 2), u 2n _ 1 = 1, and «i JV _ 1 = 0. We have to note, however, that the 
antisymmetry requirement on the function flA4j) with the spin-isospin parts being included 
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is satisfied by redefining A, u% and U2 appropriately [16(, and thus the above-mentioned 
simplicity in A, u% and U2 is destroyed when implementing the permutational symmetry 

The SA is obtained through the overlap integral of Eqs. (1A3|) and (1A4I) . The integration 
extends over all the elements of X. Because of the function 5(xn-\ — r), the coordinate 
Xjsf-i contained in Eq. (IA4h may be replaced by r. This replacement can be done as follows. 
By decomposing the matrix A into 

A -{^ a)' (A5) 
the exponent of Eq. (IA4I) is reduced to the form 

XAX = xA w x + 2^)x -r + ar 2 . (A6) 

Here the matrix A^ is symmetric and has dimension N — 2, the column vector a« is also 
of dimension N — 2, and a is equal to A^-in-i- The vector of the angular part is reduced, 
e.g., to 

U\X = Ui'x + u lN _ 1 r, (A7) 

where is a column vector which consists of the first N — 2 elements of «i. 
Now the overlap integral for the SA reads 

($LMi,(r)\F( LlLa ) L M L (ui, u 2 , A, X)} 

= ^(L'M' L £m\LM L ) [ K/ m (f) e^/exp (~xA'x ) [y L3 (u 3 x)y L4 (u 4 x)} LIML 

M' L m J 

exp ^-i^A (1) a; - a^x • r^j [y^u^x + ui N _ 1 r)y L2 (u 2 1) x + u 2N _ 1 r)] LMl ^dr .(A8) 



x 



With the change of variables from x to t, £C=t+u;r, where 

u = -B~ l a {1 \ B = A W +A', (A9) 

the overlap integral becomes 

(®lm l (r)\F (LlL2)L M L (Ui, U2, A, X)) 

^ (L' M' L £m\LM L ) J ' Y?Jr) e -^ a+a(1) "> 2 M(r)dr, (A10) 



M' L m 



where 



A<(r) = (e-^ A,t [y L3 (u 3 t + b 3 r)y Li (u 4 t + hr)] L 



M 



L 1 



with 



x e-* tAWt [y L A<h + b ir )y L2 {u ( 2 ] t + b 2 r)] LML ), (All) 
h = Ui'u + ui N _ 1 , b 2 =u^ ] uj + u 2n _ 1 , b 3 = u 3 uj, b 4 = u 4 uj. (A12) 
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To obtain the SA, we have to carry out the integration over r in Eq. (lAlOj) . The depen- 
dence of M.{r) on r is extracted using a decomposition of the type of 



y Ll (u?t + hr) = D^y Ll _ £l (u?t)y ei (b ir ) 



£1=0 



with 



D 



4vr(2L + 1)! 



(21 + 1)\(2L - 2£ + iy: 



(A13) 



(A14) 



and recoupling the angular momenta which come from the four vectors fe^r. The bra-ket 
functions depending on t lead to the overlap matrix element of the N — 2-particle system. 
After these algebraic manipulations, we get the desired formula 

(^lm l (r)\F [LlL2)LML (tlx, « 2 , A, X)) 



Li-ti h L x 
L 2 -£ 2 £ 2 L 2 

a e 12 l 



L3-4 h L 3 
A £, 4 V 



X (%3-^3i4-^4)AM A (M3,M4,^^«)|i ^ (L 1 -^ 1 L 2 -fe)AM A (M^ 1) ,M2 2 ^^ (1 ^ ^C ))> ( A15 ) 

where C is a coefficient which couples two spherical harmonics with the same argument 



C(l\l 2 ] l\ 2 ) 



f (2Zq + 1)(2Z 2 + 1) 
4vr(2/ 12 + 1) 



(ZiOZ 2 0|Z 12 0), 



(A16) 



U is a unitary Racah coefficient, and the square bracket [• • • ] stands for a unitary 9-j 
coefficient The overlap matrix element in the last line of Eq. (1A15|) can be calculated 



using Eq. (B.10) of Ref. [16 
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